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ABSTRACT 

Due to dynamical friction stellar mass black holes and neutron stars are expected to 
form high density cusps in the inner parsec of our Galaxy. These compact remnants, 
expected to number around 20000, may be accreting cold dense gas present there, and 
give rise to potentially observable X-ray emission. Here we build a simple but detailed 
time-dependent model of such emission. We find that at least several X-ray sources 
of this nature should be detectable with Chandra at any one time. Turning this issue 
around, we also ask a question of what current observational constraints might be 
telling us about the total number of compact remnants. A cusp with ^ 40 thousand 
black holes over-predicts the number of discrete sources and the total X-ray luminosity 
of the inner parsec, and is hence ruled out. Future observations of the distribution and 
orbits of the cold ionised gas in the inner parsec of Sgr A* will put tighter constraints 
on the cusp of compact remnants. 
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1 INTRODUCTION 



Theoretical calculations predict a cusp of around ~ 20000 
stellar mass black holes in the central parsec of our 
Galaxy (jMorris. 199dMiralda-Escude fc Gould. 200G| ) 

and a similar number of neutr on stars 
l|Freitag et al.. 200dHopman fc Alexander. 200^ . X-ray 
observations reveal a highly sisn ificant overabunda nce 
of transients in the same region l|Muno et al.. 2005h as 
compared with the region of ^ several tens of parsecs from 
Sgr A*, the Msmbh ~ 4 x 10^ sup er-massive black hole 
JSchodel et al.. 200jGhez et al.. 20031) in the Galactic Cen- 
tre (GC). lNavakshin fc Sunvaev (20061 ) proposed that these 
compact remnants may be accreting gas at relatively high 
rates when they happen to travel through a dense ionise d 
gas observed to exist in the GC (also see iMorris. 19931 ). 
They calculated a simple time-averaged model for X-ray 
emission from such a cusp, and concluded that the total 
emission of the cusp could be as high as 10^^ erg s~^, 
i.e. very significant observationally. 

On the other hand, a time-independent treatment does 
not do justice to the complexity of the problem. Despite 
the high total number of black holes, due to a small volume 
filling fraction of cold gas in the GC, only a few of the black 
holes will be moving within the gas clouds and possess a 
small enough relative v elocity to be visible to Chandr a. Here 
we extend the model of lNavakshin &: Sunvaev (20061 ) in two 
ways. Firstly, we allow time-dependency in the problem by 
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explicitly following realistic Keplerian orbits of the compact 
objects. Secondly, we model the formation and evolution 
of small scale accretion discs around the accretors, as such 
discs will surely form due to an excess angular momentum in 
the accreting gas. Non-circular gas orbits are also considered 
with a simplified approach. 

Our models show a large intrinsic time-dependence of 
the accretion on the compact objects and the X-ray emis- 
sion it produces. Despite that, and despite internal uncer- 
tainties of the model (exact gas orbits, circularisation ra- 
dius parameter, radiative efficiency, etc.), we feel that cer- 
tain rather robust conclusions can be drawn. In particular, 
with the ~ 20000 compact remnants expected in the cen- 
tral parsec, at least several X-ray sources with X-ray lumi- 
nosity greater than 10"^"^ erg s~^ should be present. Such 
sources, which can be called "fake X-ray binaries" for obvi- 
ous reason s, can potentially contribute to the sources ob- 
served by iMuno et al. (20051 ) in the central parsec. Con- 
versely, it appears that a cusp significantly more populous, 
i.e., with 40000 compact remnants, would over-produce the 
X-ray emission as compared to the observations, and should 
thus be ruled out. 



2 NUMERICAL APPROACH 

The basis of our setup is explained in 
iNavakshin fc Sunvaev (20061) . and hence we shall only 
briefiy discuss the main points, sparing detail for the 
improvements introduced here. 
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We assume that a stellar mass black hole travelling 
through a gas cloud or a disc with density p is captur- 
ing gas in a small s cale disc about it (see below) at the 
iBondi fc Hovle (1M4 ) accretion rate: 



capt 



: 4np 



(GMb 



(1) 



where Cs and Ad are the gas sound speed and the relative 
velocity between the black hole and the gas, respectively. 
The above picture is complicated by the presence of the 
SMBH. The area of influence of the stellar mass black holes 
is limited by the HiU's radius, m = i?(Mbh/3Msmbh)^''^, 
where R is the distance between the stellar mass black hole 
and the SMBH. This imposes a limit on the capture rate, 
given by the Hill accretion rate, Mh = 4iTr^pcB. Hence, 



Mcapt = mm [Mcapt, A/h] 



(2) 



where Mcapt on the right hand side of the equation is defined 
in eq.l. 

2.1 Time-dependent disc accretion 

The captured gas may have a net angular momentum result- 
ing in the formation of a disk around the stellar mass black 
hole. The disc size is of order of the circularisation radius for 
the gas flow, Vc, which is unknown a priori. The maximum 
value of rc is the capture radius, r-capt, which is 

GMbh 
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We thus parametrise the circularisation radius as 
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where is a parameter less than unity. The viscous time 
scale in such a disk around the stellar mass black hole is 



^visc 



1 / rcapt y 
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where a is the viscosity parameter, Sid = \/ GMhh/r'c is the 
angular velocity of the disk and h is the height scale of the 
disk. Numerically, 



tvisc = 1.5 X 10^ year ctQ^^inr^^^^ 3 



(6) 



where Td,3 is the disk temperature in units of 10 K, 
the viscosity parameter is a = O.Olao.oi, Pd is the mean 
molecular mass in units of Hydrogen mass and rc,i2 is the 
circularisation radius in lO^'^ cm. This is to be compared 
with period (the orbital time) about the SMBH: 



P = 



27r 



2900 year Rl'1M^ 



-1/2 



(7) 



Here, Q.k is the Keplerian angular frequency for the black 
hole orbiting the super-massive one: S7k = GMsmhh/ ■ 
The respective Keplerian velocity is vk ~ Ri^K- 

Thus, the gas captured in the small scale disk accretes 
on the black hole after a delay of a fraction of to a few (black 
hole around the SMBH) orbital times. The evolution of the 
disk mass is given by the rate at which the mass is added, 
Mcapt, minus the mass accreted onto the black hole, Ma.cc'- 



dt 
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(see §5.2 in iFrank et al. ('20021 ')'). Finally, the luminosity 
of the accretion flow is mod elled in the same way as 
in iNavakshin fc Sunvaev (20061 ). Namely, we write ix ~ 
eMa,ccC? , where e is given by 



0.01- 



M„., 



Mo + Ma, 



(10) 



where Mo = 0.01 is the critical accretion rate where the 
switch from th e radiatively effic ient to radiatively inefficient 
regime occurs sm et al.. 19971 ). Here we also assumed that 
X-ray emission visible in the Chandra band constitutes 10% 
of the bolometric efficiency, which would be a lower limit for 
typical spectra of X-ray binaries in their hard state. 



2.2 Orbital evolution of accretors 

We model the velocity and space dis tribution of stellar mass 
black holes as a cusp that follows the lBahcall fc Wolf ri976h 
distribution for heavier species in a mass-segregated cusp. 
This distribution results in the black hole number density 
and velocity distribution obeying power laws of the form 
and /o respectively, where. 



/o oc R 



-1/4 
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(11) 



The black hole accretion rate is calculated as 



Both space and velocity distribution s are isotropic. The 
most recent Monte-Carlo simulations iFreitag et al. (2006h 
broadly supp ort these classical results. A recent Fokker- 
Planck study l|Hopman fc Alexander. 20061 ) predicts a some- 
what steeper power-law density dependence for the black 
hole cusp, pbh oc R~^ , but we leave this level of detail for 
future investigations. 

Our black hole cusp is sha rply cut a.t R = 0.7 pc (see 
iMiralda-Escude fc Gould. 2000l ). Clearly, the artificial cut 
of the black hole distribution at the outer cusp radius is a 
crude approximation to the more c omplicated broken p ower- 
law structure of the cusp found bv lFreitag et al. (20061 ). It is 
also not entirely self-consistent as black holes do follow their 
orbits hence changing their radial position. To determine the 
significance of this we examined the structure of our model 
cusp in lO** years after it was set up. The density profile 
did in fact change. Some of black holes on eccentric orbits 
with large semi-major axes were found at radii much larger 
than the outer cusp radius assumed, as expected. However, 
the maximum change in the black hole density profile was 
no more than 30%. Furthermore, the black holes on very 
eccentric orbits will be those that accrete gas at a low rate 
unless gas in the inner parsec moves on similar eccentric 
orbits. Therefore emission from these black holes might be 
neglected in any event. 

Generating a series of orbits consistent with this space- 
velocity distribution, we randomly set the initial phases of 
the black holes along their orbits, and then we follow their 
spatial motion. We also track the instantaneous gas capture 
rate for each black hole. When one of these orbits intersects 
the disc of the Minispiral, the black hole in question starts 
capturing gas and builds up a disc around it as described in 
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Table 1. Orbital parameters of individual black hole orbits (see 
^3.11 1. The inclination of the orbit is with respect to the midplane 
of the Minispiral. The last column shows the time-averaged X-ray 
luminosity of the source. 



Black 
hole 


Inclination 

n 


Eccentricity 


Semi-major 
axis 
(pc) 


< ix > 
erg s"-"^ 


Tl 


6 


0.1 


0.1 


1.78 X 10^^ 


T2 


11 


0.2 


0.1 


2.97 X 10^'' 


T3 


1 


0.6 


0.3 


1.34 X lO^'' 


T4 


29 


0.5 


0.08 


1.28 X 10^2 



2.3 



The model for the Minispiral 

iPaumard et al. f2004l ') suggests that the Minispiral is a 
dynamical feature in a state of almost a free fall onto 
Sgr A*. We feel that the Minispiral's orbit is not likely 
to be so radial, as one would then expect Sgr A* itself 
to be accreting from the Minispiral. This would result 
in accretion rates f ar above that from the stellar winds 
l|Cuadra et al.. 2006f). and would con tradict the X-ray ob- 
servations (e.g., Baganoff et al.. 20031 ) . More realistically the 
gas in the Minispiral follows an eccentric orbit which does 
not enter the inner arcsecond (~ 0.03 parsec) of the GC. 

In our simple model, the Minispiral is modelled as half 
of a disc in a local Keplerian circular rotation around Sgr A* 
with the tot al gas mass of Mdisc = 50Mq, in accord with 
estimates in iPaumard et al. (20041 '). It extends from a ra- 
dius of 0.1 pc from the SMBH to a radius of 0.5 pc. The 
disc heightscale, H, is assumed to have a fixed ratio to 
the radius, R: H/R = 0.1. The gas density is given by 
p(R) = Mdi,,/{nR^H). 

The dynamical age of the Minispiral is a few thousand 
years. Therefore, we ran our calculations for 3000 years with 
these assumptions, and then we "remove" the Minispiral in- 
stantaneously. This is done as a rough model of time evo- 
lution of the system in the case the gas apocenter is larger 
than 0.5 parsec, so that the Minispiral would leave the inner 
0.5 parsec after a dynamical time. In >j3.4l we vary some of 
the above assumptions about the structure of the Minispiral 
to estimate sensitivity of our results to these assumptions. In 
a future work a more complicated, but unavoidably model- 
dependent dynamic of the Minispiral should be included. 



3 RESULTS 

3.1 Emission from individual black holes 

To motivate the study of X-ray emission from black hole 
cusps in this paper, we first examine the emission from in- 
dividual black holes. For simplicity of discussion in this sec- 
tion only, the "half-disk" described in H2.3I is replaced with 
a full disk, with other parameters unchanged, the only ex- 
ception being the mass of the Minispiral, which was dou- 
bled. The black holes follow Keplerian orbits ( tj2.2p that are 
characterised by the values of the semi-major axis and the 
eccentricity (see Table 1). The inclination of the orbit to 
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Figure 1. X-ray light curves of the four individual black holes 
described in i|3.1l the solid, dotted, dashed and dotted-dashed 
lines correspond to Tl, T2, T3, and T4 respectively. 



the midplane of the disc, i, is also essential in determining 
the accretion history of the black hole in question. The or- 
bital parameters of the test cases are summarised in Table 
1. Figure [T] shows the resulting light curves for the four tests 
explored. The circularisation parameter is fixed at ^ = 0.1 
for all of the tests. 

We start by looking at the most luminous case, Tl in 
Table 1, shown with the solid curve in Figure [T] Low incli- 
nation and eccentricity of the orbit ensure that the black 
hole spends all of its time inside the disc. In addition, the 
relative velocity between the black hole and the gas is small 
as the orbit of the former is close to circular. The capture 
rate is then large, resulting in a relatively high gas capture 
rate and X-ray luminosity. The few kinks in the lightcurve 
are caused by periodic variations in the relative velocity due 
to the eccentricity of the orbit. One could analyse these in 
terms of the epicyclic motion approximation for this nearly 
circular orbit. 

The black hole in test T2 is on a slightly more eccentric 
and more inclined orbit. The relative black hole-gas velocity 
is larger than in test Tl, and hence the gas capture rate (see 
eq. 1) is reduced. The black hole spends a significant amount 
of time inside the disc, but there are periods of time when it 
exits the disc through its one or the other face. Half of the 
dips in the lightcurve correspond to time spent outside the 
disc, and the other to the time when the relative velocity 
Aw reaches the highest value along the trajectory. Due to a 
relatively large size of the accretion disc (the "small scale" 
one discussed in ^2.1|l that builds up around the black hole, 
the viscous time is comparable to the duration of time spent 
outside of the disc, and hence the dips are relatively minor. 
The X-ray light curve of the source reaches a quasi-steady 
state with the luminosity Lx ~ fewx 10^* erg s~^ after ~ 10* 
years. 

In test T3 we consider a black hole on a more eccentric 
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orbit, with eccentricity e = 0.6, and a larger semi-major axis. 
The high eccentricity of the orbit results in a high relative 
velocity which limits the gas capture rate. For most of the 
orbit, Lx < 10^^ erg s~^. The luminosity of the black hole 
increases dramatically at t ~ 8000 years, when the black hole 
is at periapsis, where the relative velocity is low enough and 
the density of the gas is high enough for the luminosity to 
approach that of case T2. 

Finally test T4 is close to the worst case scenario as far 
as the gas capture rate is concerned. A high inclination and 
eccentricity orbit imply that the black hole spends little time 
inside the disc. Accumulation of gas in the small scale ac- 
cretion disc happens in a burst-like manner when the black 
hole is inside the disc. Also note that since the relative veloc- 
ity is high, gas capture radius (equation ^ is smaller than 
it is in tests Tl and T2, and hence the disc viscous time 
(equation |6} is shorter, resulting in shorter decay times for 
the "bursts" in the lightcurve. The X-ray luminosity is never 
larger than few xlO'^'^ erg s~^. 

These simple tests indicate that we may expect the 
black holes to produce a detectable X-ray emission in one of 
the two ways: (i) a few black holes may be on orbits essen- 
tially co-moving with the gas, producing a few bright point 
sources; (ii) the dim majority of high inclination and/or 
high eccentricity orbits may not produce individually bright 
sources but may be collectively bright, producing an unre- 
solvable "diffuse" X-ray emission. 

3.2 Representative cases 

Having considered the individual accretors case in the pre- 
vious section we move on to the problem of the total black 
hole cusp emission with the Minispiral model as described 
in H2.3I The upper panels in Figure 2 display the total X- 
ray luminosity of a cluster of 5000 black holes as a function 
of time for two values of the circularisation radius parame- 
ter, 1^, 0.1 and 0.001, left and right, respectively. The lower 
panels show the number of X-ray sources with luminosity 
higher than 10^^ erg s~^ for the tests shown in the panels 
just above the respective lower panels. Such sources could 
be observed by Chandra. Several conclusions can be made. 
With a larger value of = 0.1, the accretion discs around 
black holes are larger, and thus viscous times are long. As 
a result, the X-ray emission varies smoothly with time, first 
increasing as the discs are built up, and then decreasing on 
~ a thousand years time scale. Thus the sources are rather 
steady in time, and are also dim. 

For the smaller value of = 10~^, viscous times in small 
scale discs are much shorter. Therefore, the X-ray emission 
from the sources vary on much shorter time scales, i.e., of 
few years to tens of years. The sources are also brighter 
as the peak accretion rates are higher - each individual 
source shines much brighter for a shorter time, of course, 
Eis compared with the larger case. Both the upper and the 
lower panels provide us with largely independent predictions 
which may be compared to X-ray observations. 

Figure 3 shows the same experiments as Figure 2 but 
for 20000 black holes. Comparison between the two different 
values of ^ shows similar trends as before. It is interesting 
to compare the Figures 2 and 3. While the results depend 
significantly on the a priori unknown value of both low and 
high tests show same tendency of a significant luminosity 
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Figure 2. X-ray light curves (top panels) and number of individ- 
ual sources (bottom panels) where Lx < 10'^'^ erg s^^, when the 
total number of black holes in the inner parsec is 5000. The left 
and right panels correspond to = 0.1 and f = 0.001 respectively. 
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Figure 3. X-ray light curves and number of individual sources 
where Lx < 10'^'' erg s~^, when the total number of black holes 
in the inner parsec is 20000. The left and right panels correspond 
to f = 0.1 and ^ = 0.001 respectively. 



increase with increase in the number of black holes. In fact, 
the luminosity increased by a larger factor than the black 
hole number did. The number of sources above the chosen 
luminosity threshold also increased. This suggests that by 
performing tests across all reasonable parameter space for 
we should be able to find the maximum allowed number of 
stellar mass black holes in the cusp. 
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3.3 Search in the A^, ( parameter space 

Following this idea, we ran a number of models for a range 
of values of ( and for the total black hole numbers of A'^ = 
5, 10, 20 and 40 thousand. During the time period modelled, 
the results vary considerably in each test. For example, it 
is possible for just one single source to dominate the X- 
ray luminosity output of the cluster. In order to reduce and 
estimate statistical noise of the results, for each of the values 
of ( and the total black hole number considered here, the 
tests were repeated three times, each time generating a new 
random black hole orbit distribution. We then calculate the 
mean value for the observables for the three runs, and we 
also find deviations from the mean values. The averaging 
is done between time 2000 < t < 3000 years as to look 
at a state that may be similar to the present state of the 
Minispiral, given its estimated dynamical time. 

A summary of the results is presented in Table 2. We 
chose to define several quantities. The most important are 
the average total X-ray luminosity of the black hole cluster 
and the number of black holes brighter than 10^^ erg s~^, 
A'x, as such sources would have been resolved by Chan- 
dra into separate point sources. Time-dependent variants of 
these quantities were plotted in Figures 2 and 3. 

In addition to these, we defined the probabilities that 
the total luminosity of the cusp exceeds 10^^ and 10^^ erg 

as the fraction of time these conditions are satisfied: 

P{Lx> 10^') = -^ f dt, (12) 

and obviously similarly for 10"^'' erg s~^, ti and tf are 2000 
and 3000 years respectively. 

The probability of the number of individual sources 
brighter than Lx > 10'^'' being larger than 3, 10 and 20 
at a given time is defined in a similar way, e.g., 

P(Nx > 10) = / dt. (13) 

These values can be compared to the number of discreet 
X-ray sources in the inner parsec as detected by Chandra. 

Observations of the inner parsec by Chandra have 
placed upper limits on the total luminosity of s ources of 
approximately 10^^ erg s~^ (|Baganoff et al.. 20031 ) and the 
number of individual X-ray sources with a luminosity greater 
than 10^^ erg s~^ of a dozen or so (F. Baganoff, private com- 
munication). With these constraints in mind we can imme- 
diately rule out the possibility that the cusp contains 40000 
black holes. For any reasonable value of the total luminos- 
ity and the number of individual sources with Lx > 10^^ erg 
s~^ are too large compared to observations. A cusp contain- 
ing 20000 black holes is not very likely but cannot be ruled 
out completely at this time. In particular, only the larger ^ 
case is acceptable for N = 20000. Even though the average 
total cusp luminosity is ~ 10^^ erg s~^ for this test, i.e., too 
large, the probability P{Lx > 10^^) is only ~ 0.66. Cusps 
with N = 5000 or 10000 black holes is well within the limits 
imposed by observations. 



3.4 Sensitivity of results to the properties of the 
Minispiral 

We used a rather simple model for the Minispiral here ( il2.3|) . 
partially because it is not yet clear what a better model for 
this gas would be. At the moment, we do not know the ori- 
gin of this gaseous feature and the pr ecise three-dimensiona l 
distribution of gas and velocity field l|Paumard et al.. 20041 ) . 

To test sensitivity of our conclusions to the properties 
of the Minispiral, we varied several of the assumptions made 
in il2.3l Table 3 summarises these tests. In particular, in one 
of the series of tests we let the Minispiral to be 3 times 
more massive, i.e., to contain 150 M© of gas, with all other 
assumptions unchanged. In this case the luminosity of the 
cusp increases significantly at a given number of black holes, 
and even the A'^ — 10* case is too luminous for the smaller 
values of the circularisation parameter ^. Hence the upper 
limit on the number of black holes is then around 10**. 

Another likely complication is that the gas may be on 
a parabolic or an eccentric trajectory rather than a circu- 
lar one, as assumed in this paper until now. For such or- 
bits, gas velocity can be both larger and smaller than the 
local Keplerian value, depending on where exactly on the 
orbit the gas is. Observationally, the Minispiral seems to be 
clos er to the pericenter of its orbit rather than an apocen- 
ter (jPaumard et al.. 20041 ). To test the significance of non- 
Keplerian orbits, we let the gas velocity to be 1.2 and 1.4 
times the local Keplerian value in the two series of tests 
presented in Table 3. Clearly this model is not geometri- 
cally self-consistent as the half-disc we use for the Minispi- 
ral should then deform in a complicated way, but we expect 
to gain some guidance to the direction of change in the re- 
sults nonetheless. The larger gas velocity should result in 
a decrease in the number of black holes travelling at a low 
relative speed through the Minispiral, which should reduce 
the average gas capture rate (equation 1). 

Table 3 shows that increasing the gas velocity to 1.2 of 
the local Keplerian value results in a marked drop in average 
luminosity of black holes, to the point that a cusp containing 
40000 black holes cannot be ruled out for the larger value of 
the circularisation parameter, ^ = 0.1. On the other hand, 
realistically, we would expect the sources to have a distribu- 
tion in values of (, and we would thus expect a fair number 
of sources to have zeta ~ 0.01 or less, which would then be 
ruled out. 

Moving on to the gas velocity = IAvk, we note a fur- 
ther drop in the luminosity of the cusp and the number of 
sources detectable by Chandra (Table 3). Even the A'^ = 
40000 is allowed. We however feel that this model strongly 
over-simplifies the situation in the Galactic Centre. The 
Minispiral should be a feature bound to Sgr A* (or else the 
fact that it is crossing the innermost region of the parsec now 
would be a coincidence), and hence it is rather implausible 
that the gas is on a radial trajectory for which one would 
have V = V^vk ~ 1.4t;K. We believe the case of i; = 1.2i;k 
is the one that represents the reality better. 

Finally, the last entries in Table 3 are the tests with 
the Minispiral mass of M = 150 Mq and v = 1.2i;k. A 
cusp of 40000 black holes is clearly inconsistent with the 
observations, whereas the A'^ = 20000 is constrained but not 
completely ruled out. 
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3.5 A neutron star cusp 

The mechanism that produces an overabundance of black 
holes in the Galactic Centre also applies to neutron stars, as 
they are also heavi er than an average star. Simulations by 
iFreitag et al. f2006l) predict that a cusp of neutron stars will 
have a number density profile quite similar to that of stel- 
lar mass black holes. As far as our investigation goes, these 
neutron stars would be accreting gas in a similar manner to 
the black holes, with modifications only due to the smaller 
mass (we set Mns ~ 1.4 Mq) and the existence of a surface. 

We approximate the emission of an accreting neutron 
star as a black body one with temperature Tns from the 
surface area where i?ns = 10 km is the radius of the 

neutron star. Most of the radiation fiux will be emitted at 
wavelengths corresponding to photon energy E = SkTns'- 

E « 0.8 Lgf KeV , (14) 

where L34 is the X-ray Luminosity in units of lO"** erg s~^. 
Due to the quite large absorbing colum n density to the GC, 
A'h ~ 10 cm~ ^ (|Baganoff et al.. 2003l \ soft X -ray emission 
below ~ 1 keV is practically unobservable. We hence set the 
minimum observable total X-ray luminosity of a neutron star 
to be lO'^* erg s~^, rather than 10''^ erg s~^ for the black 
hole case. 

We assume that the radiative efficiency of an accreting 
neutron star with negligible magnetic fields is constant at 
e ~ 0.1, in which case the luminosity of an accreting neutron 
star is simply 



0.1 Afa, 



(15) 



Note that if magnetic fields are non-negligible, then the 
"propeller effect" may reduce the X-ray luminosity of neu- 
tron stars (|Menou et al.. 19991 ). but we leave these compli- 
cations for a future more extended study. 

With these modifications, we can use the machinery 
developed in ij2] to calculate the X-ray emission from the 
cusp containing neutron stars. The results are presented in 
a way identical to the black hole cusp in Table 4. As with 
the black hole cusp case, the A'' = 40000 case is strongly 
ruled out on account of too large a number of detectable 
point sources, A'x, and the total X-ray luminosity of the 
cusp. The A'' = 20000 cusp also appears to be a bit too 
high in terms of both the total luminosity and the number 
of sources. The data may still tolerable this, especially if 
one takes into accou nt the possibility of the propeller effect 
jMenou et al.. 19991 ) neglected here. 



4 CONCLUSIONS 

Stellar mass black holes and 
dieted to clutter the central 



neutron stars are pre- 
parsec of our Galaxy 



elaborated study, where a time-dependent disc accretion 
onto the compact sources was considered. Also, we took 
into account the fact that at low accretion rates, the 
radiative effic iency of black ho les appears to be drastically 
reduced ('e.g.. lEsin et al.. 1997I ). and we used a Monte-Carlo 
like approach to randomly initialise the cusp of compact 
remnants. 

The main effort in our paper was to set the upper limit 
on the number of compact remnants. Whereas the models 
have internal uncertainties, such as the value of circularisa- 
tion parameter ^, and observational uncertainties (the Mass 
and precise orbit of the Minispiral) , a cusp of black holes or 
neutron stars with A'^ ^ 40, 000 seems to be strongly ruled 
out. A cus p with ^ 20, OOP black holes, as theoreticall y 
predicted (|Freitag et al.. 200dHopman fc Alexander. 20061 ). 
is broadly consistent with the data. Future efforts should 
improve these upper limits. 

On the basis of our calculations, we think it is quite 
realistic tha t some of the X-ray sources visible in th e central 
parsec (e.g.. lBaganoff et al.. 200jMuno et al.. 20051 ) may be 
isolated black holes and neutron stars accreting gas from the 
Minispiral. Such sources should be preferentially found close 
to the Minispiral if viscous time is short {( is small). In ad- 
dition, binary systems containing a black hole and a normal 
low mass star can also accrete gas in roughly the same way 
as we calculated here. In the case of low values of circulari- 
sation parameter, the size of the disc around the primary 
(the black hole) can be smaller than the size of the binary 
itself. Thus, these systems may appear as "fake X-ray bina- 
ries" , where the gas supply comes from outside rather than 
from the low mass secondary. Observational signatures of 
such systems might be warped and out of binary plane ac- 
cretion discs, "too short" or "too weak" accretion outbursts 
for the size of the binary. 

We acknowledge stimulating discussions with Mark Fre- 
itag, Fred Baganoff and Rashid Sunyaev. 
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Table 2. Characteristics of black hole cusp averaged between 2000-3000 years (see i|3.3| l. 



Number" 
[103] 




[lO^'^erg s^i] 


P(Lx > W^^Y 


P(Lx > 10^^)'= 


< Nx >f 


P(Afx >20)9 


P(Afx >10)'* 


P(Afx >3)' 


5 


0.001 


1.85 ± 0.51 





0.53 ± 0.5 


2.76 ± 0.75 





0.00 ± 0.00 


0.33 ± 0.8 


5 


0.01 


2.75 ± 1.33 


0.42 ± 0.2 


0.72 ± 0.1 


4.2 ± 0.4 








0.63 ± 0.13 


5 


0.1 


0.73 ± 0.6 





0.24 ± 0.3 


2.23 ± 0.2 








0.18 ± 0.11 


10 


0.001 


15.91 ± 2.70 


0.33 ± 0.13 


0.90 ± 0.07 


4.44 ± 0.28 








0.71 ± 0.06 


10 


0.01 


1.03 ± 0.24 





0.45 ± 0.16 


5.36 ± 0.26 








0.90 ± 0.04 


10 


0.1 


0.51 ± 0.28 





0.31 ± 0.22 


4.53 ± 1.52 





0.01 ± 0.00 


0.45 ± 0.21 


20 


0.001 


19.97 ± 5.87 


0.57 ± 0.05 


0.99 ± 0.01 


8.87 ± 0.81 





0.19 ± 0.11 


1.00 ± 0.00 


20 


0.01 


9.67 ± 1.94 


0.34 ± 0.14 


0.99 ± 0.00 


12.26 ± 1.40 


0.01 ± 0.01 


0.69 ± 0.17 


1.00 ± 0.00 


20 


0.1 


9.88 ± 4.49 


0.34 ± 0.22 


0.66 ± 0.23 


13.11 ± 1.25 


0.01 ± 0.00 


0.74 ± 0.16 


1.00 ± 0.00 


40 


0.001 


40.37 ±1.75 


1.00 ± 0.00 


1.00 ± 0.00 


19.40 ± 0.37 


0.34 ± 0.05 


1.00 ± 0.00 


1.00 ± 0.00 


40 


0.01 


100.53 ± 62.13 


0.77 ± 0.08 


0.99 ± 0.00 


24.66 ± 0.50 


0.85 ± 0.05 


1.00 ± 0.00 


1.00 ± 0.00 


40 


0.1 


6.50 ±1.83 


0.20 ± 0.14 


1.00 ± 0.00 


24.59 ± 2.18 


0.72 ± 0.14 


1.00 ± 0.00 


1.00 ± 0.00 



The coUumns list: 

Total number of black holes in the cusp 
^ Circularisation parameter ( ^2.11 1 

Time-averaged luminosity of the cusp 
d-e Probability that the total luminosity of the cusp is greater than 10^® or 10"^^ erg s~^, respectively. 
^ Average number of sources with X-ray luminosity greater than lO^^erg s~^{Nx) 
Probability that A^x is greater than 20, 10 and 3, respectively 

Muno M. P., Pfahl E., Baganoff F. K., et al., 2005, ApJL, 622, 
L113 

Nayakshin S., Sunyaev R., 2006, ArXiv Astrophysics e-prints 
Paumard T., Maillard J.-P., Morris M., 2004, A&A, 426, 81 
Schodel R., Ott T., Genzel R., et al., 2002, Nature, 419, 694 
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Table 3. Same as Table 2, but for different models of the Minlspiral 



Number 

[103] 


c 


<Lx > 
[lO^Sg^g g-ij 


P(Lx > 10^^) 


P{Lx > 1035) 


< ATx > 


P{Ny, >20) 


P(Arx >10) 


P(iVx >3) 


3 Mgas'^ 
5 


0.01 


28.13 ± 11.54 


0.62 ± 0.15 


0.84 ± 0.08 


5.81 ± 0.95 





0.03 ± 0.02 


0.84 ± 0.06 


5 


0.1 


1.53 ± 4.45 





0.67 ± 0.40 


7.04 ± 0.47 





0.06 ± 0.04 


1.00 ± 0.00 


10 


0.01 


35.21 ± 12.01 


0.63 ± 0.17 


1.00 ± 0.00 


9.24 ± 0.93 





0.36 ± 0.13 


1.00 ± 0.00 


10 


0.1 


2.06 ± 0.74 





0.56 ± 0.16 


14.12 ± 0.51 


0.02 ± 0.00 


0.98 ± 0.01 


1.00 ± 0.00 


20 


0.01 


24.53 ± 6.02 


0.77 ± 0.13 


1.00 ± 0.00 


18.95 ± 1.25 


0.40 ± 0.16 


1.00 ± 0.00 


1.00 ± 0.00 


20 


0.1 


5.47 ± 0.41 


0.01 ± 0.00 


1.00 ± 0.00 


29.38 ± 2.06 


0.90 ± 0.07 


1.00 ± 0.00 


1.00 ± 0.00 


1.24 

5 


0.01 


93 + 55 


02 + 01 


1 4 + 07 


1 85 + 36 








1 + 04 




0.1 


0.05 ± O.O.'i 








0. 11 ± 0.20 











10 


0.01 


0.45 ± 0.10 





0.13 ± 0.05 


3.47 ± 0.42 








0.49 ± 0.13 


10 


0.1 


0.38 ± 0.25 





0.17 ± 0.12 


2.45 ± 1.42 








0.31 ± 0.22 


20 


0.01 


4.86 ± 1.63 


0.14 ± 0.05 


0.43 ± 0.15 


6.65 ± 0.08 





0.03 ± 0.01 


0.91 ± 0.04 


20 


0.1 


0.35 ± 0.07 








4.24 ± 0.08 








0.76 ± 0.07 


40 


0.01 


40.65 ± 1.66 


0.56 ± 0.14 


0.94 ± 0.04 


13.42 ± 0.23 


0.02 ± 0.01 


0.74 ± 0.06 


1.00 ± 0.00 


40 


0.1 


1.05 ± 0.37 





0.18 ± 0.13 


10.27 ± 0.62 





0.44 ± 0.80 


1.00 ± 0.00 


1.4 
5 


0.01 


n 07 + n 02 








n 71 + n 1 c; 











5 


0.1 


0.02 ± 0.00 








0.37 ± 0.26 











10 


0.01 


0.12 ± 0.04 








1.17 ± 0.21 











10 


0.1 


0.05 ± 0.01 








0.83 ± 0.32 











20 


0.01 


0.32 ± 0.48 





0.03 ± 0.02 


3.66 ± 0.43 








0.49 ± 0.13 


20 


0.1 


0.03 ± 0.00 








0.11 ± 0.05 











40 


0.01 


0.64 ± 0.14 





0.18 ± 0.09 


5.37 ± 0.60 





0.04 ± 0.01 


0.72 ± 0.07 


40 


0.1 


0.12 ± 0.03 








1.34 ± 0.36 








0.01 ± 0.01 


3 Mgas & 
20 


1.24 
0.01 


9.79 ± 0.93 


0.34 ± 0.06 


1.00 ± 0.00 


15.47 ± 0.25 


0.07 ± 0.02 


0.98 ± 0.01 


1.00 ± 0.00 


20 


0.1 


2.89 ± 1.03 


0.09 ± 0.06 


0.86 ± 0.10 


13.72 ± 0.25 


0.04 ± 0.01 


0.91 ± 0.03 


1.00 ± 0.00 


40 


0.01 


12.17 ± 0.96 


0.59 ± 0.02 


1.00 ± 0.00 


25.02 ± 0.78 


0.90 ± 0.05 


1.00 ± 0.00 


1.00 ± 0.00 


40 


0.1 


6.43 ± 1.76 


0.29 ± 0.10 


0.93 ± 0.05 


19.63 ± 1.82 


0.50 ± 0.16 


1.00 ± 0.00 


1.00 ± 0.00 



" Mass of the Minispiral has been tripled to 150 

^ Velocity of the gas has been increased to 1.2 vj<^, where is the local keplerian velocity 
Velocity of the gas has been increased to 1.4 vk 

Mass of the Minispiral has been tripled to 150 Mq and the velocity of the gas has been increased to 1.2 vj^ 
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Table 4. Same as Table 2 but for a neutron star cusp. 



Number" 
[10«] 




[lO^^crg s"^] 


P(Lx > lO^^)"* 


P(Lx > 1035)e 


< A^x >^ 


P(Afx >20)9 


P(Arx >10)'^ 


P(Afx >3)" 


5 


0.001 


3.23 ± 5.24 


0.08 ± 0.06 


0.79 ± 0.11 


2.84 ± 0.43 








0.32 ± 0.10 


5 


0.01 


2.43 ± 0.76 





0.67 ± 0.12 


2.78 ± 0.81 








0.38 ± 0.16 


5 


0.1 


1.44 ± 0.18 





0.85 ± 0.11 


2.90 ± 0.67 








0.34 ± 0.20 


10 


0.001 


4.17 ± 0.90 


0.02 ± 0.01 


0.99 ± 0.01 


6.25 ± 0.45 





0.03 ± 0.02 


0.94 ± 0.02 


10 


0.01 


5.89 ± 1.38 


0.18 ± 0.08 


0.99 ± 0.01 


4.53 ± 0.34 








0.72 ± 0.10 


10 


0.1 


6.69 ± 3.59 


0.30 ± 0.21 


0.83 ± 0.12 


3.32 ± 0.15 








0.47 ± 0.05 


20 


0.001 


8.13 ± 1.34 


0.27 ± 0.11 


1.00 ± 0.00 


11.52 ± 1.15 





0.57 ± 0.15 


1.00 ± 0.00 


20 


0.01 


9.73 ± 1.73 


0.32 ± 0.12 


1.00 ± 0.00 


7.32 ± 0.30 





0.11 ± 0.08 


0.97 ± 0.02 


20 


0.1 


7.17 ± 2.15 


0.27 ± 0.19 


1.00 ± 0.00 


9.00 ± 0.35 





0.19 ± 0.05 


1.00 ± 0.00 


40 


0.001 


41.60 ± 12.71 


0.99 ± 0.01 


1.00 ± 0.00 


21.99 ± 0.52 


0.69 ± 0.07 


1.00 ± 0.00 


1.00 ± 0.00 


40 


0.01 


18.62 ± 4.46 


0.87 ± 0.07 


1.00 ± 0.00 


15.82 ± 0.48 


0.04 ± 0.02 


0.99 ± 0.00 


1.00 ± 0.00 


40 


0.1 


17.91 ± 3.74 


0.86 ± 0.10 


1.00 ± 0.00 


21.34 ± 0.34 


0.66 ± 0.05 


1.00 ± 0.00 


1.00 ± 0.00 



<^ — see Table 2 caption. 

f Average number of neutron stars with Lx greater than lO'^'^erg s~^{Ny^). Only sources with Lx > 10^* are visible, see i|3.5l 
5 ~ ' see Table 2 caption. 



